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We study the potential of tlie LHC to discover the signal of Dark Matter associated production 

Qh| with a photon induced by a dimension six effective operator, including NLO QCD corrections. We 

Q\ • investigate the main backgrounds from SM, i.e. Z boson and a photon associated production with 
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o 



s 



invisible decay of Z boson, and Z boson and a jet production with the jet misidentified as a photon. 
We find that the pj^ distributions of the backgrounds decrease faster than that of the signal with 
increasing of the transverse momentum of the photon. The rj'^ distributions of the backgrounds are 



%^j I We find that the pj, distributions of the backgrounds decrease faster than that of the signal with 

OO 

o 

^\i ■ almost flat in the full range of rj'^ . In contrast, the signal lies mainly in the central region of r]'^ . 



These characteristics may help to select the events in experiments. We show that in the parameter 
space allowed by the relic abundance constraint, which we have calculated at the NLO QCD level, 
the LHC with v5 = 7 TeV may discover this signal at the 5a level after collecting an integrated 
'^ . luminosity of 1 fb . On the other hand, if this signal is not observed at the LHC, we can set a 

lower limit on the new physics scale at the 3a level. 
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I. INTRODUCTION 

Astrophysical and cosmological observations have confirmed the existence of Dark Mat- 
ter (DM) in our universe [0] and the density of DM is much larger than that of the visi- 
ble matter 0. The relic abundance of DM favours a weakly interacting massive particles 
(WIMP). This kind of DM has been extensively studied in the literatures. Since there are 
no candidates of this kind of DM in the Standard Model (SM), any discovery of the signal 
of DM imply new physics. 

In the region of DM accumulating, DM can annihilate into SM particles, such as photons, 
electrons and positrons. These produced particles can propagate through the interstellar 
space and be detected by experiments on the earth, such as PAMELA[||, ATICIQ, HESSp 
and Fermi LATp. The detection of these signals is not conclusive evidence for DM, since it 
depends on the assumptions of the distribution of DM and the propagator model. Besides, 
other astrophysical interpretation can not be excluded. Another way to look for the signal of 
DM is to measure the recoil energy of nuclei caused by the elastic scattering of a WIMP off a 
nucleon, such as the experiments of DAMA[0], CDMS||, CoGeNT|| and XENON^Ol. The 



DAMA and CoGeNT experiments favour a light DM with a mass around 10 GeV. The CDMS 
and XENON experiments set upper limits on the WIMP and nucleon spin-dependent and 
spin-independent cross sections if the mass of the WIMP ranges from 10 GeV to 1000 GeV. 
These experiments are passive and much time was spent in waiting for the collision with 
the DM. More active approach is to produce the DM in the laboratory directly, such as 
the Large Hadron Collider (LHC) if DM exists and has interactions with the SM particles. 
There are a lot of studies to search for DM at the LHC in varies of DM models ||TT|-|28|. 

Because the LHC is a proton-proton collider, the QCD correction should be considered 
for any process if people want to make a reliable prediction. In this work, using model 
independent method, we investigate the possibility of discovering the DM in associated 
production with a photon induced by a dimension six effective operator at the next-to- 
leading (NLO) order QCD level, since this signal is clear and suffer from little backgrounds 
from the SM. 

This paper is organized as follows. In section |I|, we describe the dimension six effective 
operator for this process. In section |T|, we calculate the relic abundance induced by this 
effective operator and find the allowed region for the mass and couplings of the DM. In 



section |V[ we present the details of the NLO QCD corrections to the associated production 
of the DM and photon and discuss the dependence of the K-factor on the mass and couphngs 
of the DM. In section |V], we calculate the backgrounds in SM and analyze the discovery 
potential at the LHC. Conclusion will be given in section |V1[ 

II. EFFECTIVE OPERATOR 

The DM studied in this work is a Dirac fermion, denoted by x- It is a singlet under the 
SM gauge group 5'f/(3)c x SU(2)l x U(1)y- As a result, this kind of DM does not participate 
in any strong or electroweak interactions. In addition, we assume it couples with the SM 
quarks in the form 

o = ^,m{xxi (1) 



where A is the new physics scale above which new particles should appear. This four-fermion 
operator has been considered in Refs. [^, p9l-|31| in discussing DM interactions. However, we 
will focus on the signal of DM and photon associated production at the LHC which has not 
been carefully investigated before. Moreover, we calculated the NLO QCD corrections for 
this process whose effects are important for research at the LHC. Note that the new physics 
scale A can be viewed as the remnant of integrating the propagator between the SM particles 
and DM. Therefore, this operator is valid only if A > -m, Vs, where m is the mass of DM 
and Aj/s the center of mass energy of the collision. Generally, it is possible that a/S > A at 
the LHC. However, the luminosity for the process drops very fast with the increasing of Vs 
at the LHC with yS = 14 (7) TeV, which can be seen in Fig.l. Thus, we relax this limit in 
practical numerical calculation, where we set the default value of A=500 GeV. 

III. RELIC ABUNDANCE 

The DM relic abundance is a precision observable in cosmology. The DM we choose 
to study can contribute to the relic abundance of cold Dark Matter (CDM). Combining 
WMAP data with the latest distance measurements from Baryon Acoustic Oscillations in 
the distribution of galaxies and Hubble constant measurements gives the constrain [R^] 



QcDMh^ = 0.1123 ± 0.0035, (2) 
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FIG. 1: Luminosity plots for uu initial states at the LHC with v5 = 14 (left) and 7 (right) TeV. 



where flcDM is the CDM energy density of the Universe normahzed by the critical density 
and h = 0.710 ± 0.025 is the scaled Hubble parameter. 

The relic abundance can be calculated from the total annihilation cross section of DM. 
First, we give the LO total annihilation cross section as 

K^ s — 4m^ 



a'^^v = N,Nf- 



(3) 



A4 Svr ' 
where v is the relative velocity between the DM. A^^ and Nj are the numbers of color and 



flavor of quarks, respectively. This result agrees with that in Ref. |29 



The NLO corrections to the total annihilation cross section comprise of two parts: the 
one-loop virtual corrections and real gluon emission corrections. The results of virtual 
corrections are 
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where D^ = {Airfi'^ /sY/ r{l — e) and o"^"''^ is the n-dimensional (n = 4 — 2e) LO total 
annihilation cross section 
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In obtaining the above results, we have renormalized the effective operator in the MS scheme. 
The results of real corrections are 
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Combining the two parts we get the NLO total annihilation cross section 
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The DM is moving at nonrelativistic velocities(t' ^ 1) when freezing out. Thus we can 
expand 

a'^«t; = a + bv^, (8) 

where 



a = 0, 
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b = K'^'-N^Nf--, (9) 



in which K"'"' is the K-factor of the DM annihilation cross section 
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The freezeout epoch Xf = m/Tf is determined by| 



Xf = In [0A56bmpim (g/gl/^)] - ^ In {in [0A56bmpim (g/gl/^)] } , (11) 

where mpi = 1.22 x 10^^ GeV is the Planck mass, g and g^ is the number of relativistic 

degrees of freedom and effective number of relativistic degrees of freedom at the freeze-out 

temperature Tj, respectively. The result for the relic abundance is 

, 1.07 X 10^ GeV~^xl , , 

fix^ ~ m -■ (12) 

^{g*s/g* )mpib 

For most of the history of the universe all particle species had a common temperature, and 
g^:S can be replaced by g^:. Requiring the DM relic abundance is in the 2a region around 
the observed central value, the new physics scale of the effective operator is determined by 
the mass of the DM. We show this relation in Fig.^. When the parameters {m, A) satisfy 
this relation, the NLO QCD K-factor is nearly 1.4. Comparing the LO result (blue band) 
and the NLO result (red band), we find that the NLO QCD correction increases the new 
physics scale by about 10%. Since we do not hold the point of view that the abundance of 
the WIMP is determined only by this one kind of DM, the regions below the red band are 
all allowed. 



IV. NLO QCD CORRECTIONS TO DM AND PHOTON ASSOCIATED PRO- 
DUCTION 

According to the operator in (|l|), the DM can be pair produced at hadron colliders. How- 
ever, because DM can not decay into SM particles, such processes give just missing energy 
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FIG. 2: Relation between the mass of DM and the new physics scale. The relic abundance is 
required be to in the 2a region around the observed central value. The blue band is the LO result. 
The red band is the NLO result. In this figure, we choose k = 1, a^ = 0.118 and Nf = 5. 

and no observable signals. As a result, we have to consider the process of DM associated 
production with a photon or a jet [pi], |30| , [3l[] . Though the cross section of associated pro- 
duction with a photon is less than that of associated production with a jet, the signal of this 
process is characteristic and bring with less backgrounds. 



A. LO calculation 

We start from the leading order (LO) calculation. The LO Feynman diagrams for this 
process 

q{pi) + q{p2) -^ xiPs) + xiPi) + l{P5) (13) 

are shown in Fig.^ 

Because the DM does not interact with SM particles except for the quark fields in the 
operator (P, their contributions to the cross section can be factorized as 



|-MDMr = 2(s34-4m2) 



(14) 
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FIG. 3: LO Feynman Diagrams, 
where Sij = {pi +Pj)^. Therefore, the born matrix element can be written as 

-2M3ti5 + Ml (-ti5 - t25) + 2M2t25 , . .,^. 

Mb = Mdm, (15) 

tl5't25 

where the Mj,i = 1,2,3 represent the standard matrix elements in our calculation which 
are defined as 

Ms = v{p2)u{pi)p'iei,{p5), 

Ms = v{p2)u{pi)pi^e^{p5). (16) 

The spin and color summed and averaged Born matrix element squared is 

\Mb\ - ^TTl ri- \-^dm\ , [i-l) 



where tij = {pi — pjY and a = e'^/Air. Then the LO partonic cross section i 



is 
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^B = ^J dTslMsl^ (18) 

in which Fs is the three particle final states phase space. After convoluting with the parton 
distribution functions (PDFs) G'q(g)(x), we obtain the LO cross section 

(Tb = j dxidx2[Gg/p{xi)Gg/p{x2) + {xi f> Xs)]^-^. (19) 



B. NLO results 

The NLO QCD corrections involve the one-loop virtual gluon effects and contributions 
of real gluon and quark or antiquark emissions to leading order processes. To deal with 
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FIG. 4: Feynman Diagrams for one-loop virtual corrections. 



ultraviolet (UV) and infrared (IR) (soft and collinear) divergences in our computation, we 
use n = 4 — 2e dimensional regularization to regulate these divergences, and all divergences 
appear as 1/e* with i = 1, 2. 

The virtual gluon corrections to DM and photon associated production consist of the 
self-energy, vertex and box diagrams, as shown in Fig.§. The UV divergences are canceled 
between the loop diagrams and counterterms which are determined in on-shell renormaliza- 
tion scheme for external fields and in MS scheme for coupling constants. The final virtual 
gluon corrections to the partonic cross section are 

1 
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where C^ = 1(1 + e)[(47r/x^)/si2]' and 
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We can also write Eq. (EO) as 
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FIG. 5: Sample Feynman Diagrams for gluon-gluon initial states contributions. 
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FIG. 6: Feynman Diagrams for a real gluon emission. 



with 



da„ 
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The first term in A^ results from the MS renormalization of the four-fermion operator. 
The IR divergences remain after renormalization and have a structure predicted by other 
methods. They will cancel the IR divergences coming from real corrections. 
A further contribution at hadron colliders arises from the process 



9 + 9^X + X + 7 



(25) 



which is illustrated in Fig.|^. The effects of these loop-induced gg diagrams are one order of 
as higher than that of qq diagrams, but they may be still important due to the large gluon 
PDFs at the LHC However, the triangle and box diagrams in Fig.^ are both vanishing 
because of the color structure and charge conjugation symmetry, respectively. 
The Feynman diagrams for the real gluon emission process 



q{pi) + q{p2) -^ xiPs) + xiPi) + 7(P5) + gipe) 



(26) 



are shown in Fig.^. When performing the final states phase space integration, one encounters 
the soft and coUinear singularities. We use the two cutoff phase space slicing method to 



separate the singular regions and perform the integration analytically in these regions [p4 . 
Therefore, the real corrections are divided into three parts, i.e., 

d&r = da^ + da^^ + daf^, (27) 

where a^ and a^'" denote the contributions from soft and hard collinear regions, respectively. 
The soft regions are defined by the energy of the emitted gluon Eg < (5sy^ii2/2, where 6s 
is the soft cutoff parameter. The collinear regions are determined according to whether 
the Mandelstam variables tjg = {Pi ~P&Yi with i = 1,2, satisfy the collinear condition 
l^iel < 6cSi2, where 6c is the collinear cutoff parameter. The hard non-coUinear part cr^^'-' is 
finite and can be computed numerically. 

The partonic cross section in soft regions can be factorized as 

da^ = {A-nasiil')daB I dS^eik, (28) 

where dS is the integration over the phase space of the soft gluon 

dS = ^^J^3_,, ^ ' '" dE,El-'^dn2-2e (29) 



and the Eikonal factor $eifc is the amplitude squared in the soft limit apart from |A^bP, 
expressed as 

^eik = Cf-^. (30) 

After integration over the soft gluon phase space, we get 

da^ = das^a (4 + — + ^o] , (31) 

where 

Al = 2Cf, Af = -ACf In 5„ A^ = Cf U In' <5, - ^"j . (32) 

In hard collinear regions of this process, the momentum of the gluon emitted from initial 
partons become collinear to the beam line. In this limit, the four-body matrix elements are 
approximated as follows: 



\Mr\'' ^ {A'Kasfit)\MB\ 



-2P,,iz,e) ^-2Pgg{z,e) 



(33) 



ZtiQ Zt26 

in which, z represents the fraction of initial partons' momentum carried by q{q). Pij{z,e) 
are the unregulated splitting functions in n-dimensions which can be related to the usual 
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Altarelli-Parisi splitting kernels as Pij{z,e) = Pij{z) + eP^j{z). In our case, 

1 + ^2 



Paa(z) = Cf 



qq 



1- z 



PJ^) 



-Cp(l-z). 



At the same limit, the four-body phase space can be written as 

(Any 



dV 



4 coll 
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167r2r(l - e) 
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Therefore, we obtain 



da, 



HC 



= das^C, ( — j 6^ '[Pgg{z, e)Gg/p{xi/z)Gg/p{x2) 
+ Pqg{z, e)Gg/p{xi)Gg/p{x2/ z) + (xi ^ X2)] — 



dxidx2- (36) 

z \ z 

To factorize the collinear singularity into the PDFs, we use scale dependent PDFs in the 
MS convention: 



Gb/p{x,fif) = Gb/p{x) + 



as r(l — e) / 47r/x, 



^ dz 
—PUz)Ga/p{x/z). (37) 

z 



2tx r(l - 2e) \ [x) 

Now, we replace Gq(^q)ip in the LO hadronic cross section (|19D and combine the result with 
the hard collinear contribution (|36D . The resulting 0{as) expression for the initial state 
collinear contribution is 



a.. 



da"" = daB—G,\Gg/p{xi,i2f)Gg/p{x2,iJ,f) + Gg/p{xi,iJ,f)Gq/p{x2,fif) 



Af (a -)■ ag) 



a=q,q 

+ {xi -^ X2) >dxidx2- 



+ A^Q^^a -^ ag) Gq/p{xi, ^if)Gq/p{x2, /i/) 



with 



(38) 



ATiq^qg) = Cf(21n5, + 3/2), 
Al^q^qg) = ATiq ^ qg) In ('-^). 

\fXf/ 



The G functions are given by 



with 



Gb/p{x,fif) = 22 I —Ga/p{x/y,fif)Pba{y) 



PUy) = PUy) In ( 5^—y-'-^] - PbM- 



y /^/ 
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(41) 
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FIG. 7: Feynman Diagrams for a quark or antiquark emission. 

A complete real correction includes also the (anti) quark emitted processes, as shown in 
Fig.|^, such as 



g{pi) + q{p2) -^ xiPs) + xiPi) + liPb) + q{p&)- 



(42) 



Their contribution can be obtained from the results of processes in ( p6D by crossing sym- 
metry. Note that, when separating the singular phase space regions for these processes, one 
only needs to deal with the coUinear divergences which can be totally absorbed into the 
redefinition of the PDFs in (^) . 

Finally, the NLO cross section for the process pp — ?■ XXI is 



a 



NLO 



coll 



dxidx2\[Gq/p{xi, ^,f)Gq/p{x2, ^if) + {xi ^ X2)]{aB + 0-^ + 0-*^ + o"f'^)| +cr' 
+ ^ / dxidx2[Gg/p{xi, i^f)Ga/p{x2, fif) + {xi O X2)](fr'-' {ga -^ XX7«), (43) 



a=q,q 



where G in dr^{ga — )■ XXI'^) suggests that the phase space integration is performed in the 
non-collinear regions. We have checked that 



Al + Al 



0, A\ + Af + 2Anq^qg) = Q. 



(44) 



Therefore there are no singularities left now. 



Numerical results 



In this subsection, we give the numerical results for the cross sections for DM and pho- 
ton associated production at the LHC In numerical calculation, we choose the CTEQ6L1 
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(CTEQ6M) PDF sets|3^ and the corresponding strong coupling as for the LO (NLO) cal- 
culations. The default factorization and renormalization scales, fif and fir, are set as 2m. 
We choose the input parameters (m, A) = (150 GeV, 500 GeV) and k = 1 unless otherwise 
specified which are allowed by the relic abundance constraint. The kinematic cuts 

pI > 100 GeV, 

\ri^\ < 2.4, 
pip''' > 100 GeV, (45) 

are applied in our numerical calculation. Here p^^^^ is the missing transverse momentum, 

defined as 

{pZ, no iets in the final states, 
II (46) 

Py , with jets in the final states, 

where p^^ is the transverse momentum of the system of the DMs. Jets are defined by the 
requirements p^^ > 20 GeV and \ri^'^^\ < 2.5. In order to avoid QED collinear divergences, 
we also require the photon to be isolated by the prescription |]36[ 



5:^"<^j(\^^)- (*^) 

where R = a/A(/)2 + Ar/^ and Rq = 0.4. 

In Fig.| we show the cutoff parameter dependence of the NLO cross sections. The 
contribution of three-body final states includes the Born cross section, one-loop virtual 
corrections, soft and collinear limits of the cross section of four-body final states. The 
contribution of four-body final states consists in the cross section of four-body final states 
with the singular regions of the phase space sliced. The change of the NLO result is very 
slow, especially in the region 6s < 10^^, which indicates that it is reasonable to use the two 
cutoff phase space slicing method. 

In Fig.^ we show the m and A dependence of the LO and NLO cross sections. The LO 
(NLO) cross sections decrease from 1097 (1082) fb to 922.3 (936.3) fb as m increases from 
130 GeV to 200 GeV. The corresponding K-factor, defined as the ratio of the NLO cross 
sections to the LO ones, varies from 0.99 to 1.02. The LO (NLO) cross sections decrease 
from 1044 (1041) fb to 65.24 (70.71) fb as A increases from 500 GeV to 1000 GeV. The 
corresponding K-factor varies from 1.00 to 1.08. The NLO QCD corrections are modest. 

13 



8000- 







- 


6000- 




4-body 

NLO 


4000- 




3-body 


■ 


2000- 




- 


0- 




■- 


2000- 




- 


4000- 




- 


6000- 




- 


8000- 





-I-I-TT 



1100 

1050 

1000 

950 




FIG. 8: Dependence of the NLO cross sections for the DM and photon associated production at 
the LHC on the soft cutoff 6s with 6c = 6s/ 50. 
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FIG. 9: Dependence of the LO and NLO cross sections for the DM and photon associated produc- 
tion at the LHC on the DM mass and the new physics scale A. Also shown is the K-factor. 

However, the dependence of the NLO cross section on the factorization scale fip and renor- 
malization scale fin is significantly reduced, as shown in Fig.|TU|. This makes the theoretical 
prediction much more reliable. 
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FIG. 10: Dependence of the LO (NLO) cross sections for the DM and photon associated production 
at the LHC on the factorization scale ^p and renormalization scale ^u/j. 

V. BACKGROUNDS AND DISCOVERY POTENTIAL 



The dominant SM backgrounds for this process include the processes qq — ?■ Z{^ vv) +7 
and qq — t- Z{-^ vv) + j with the jet misidentified as a photon. The NLO QCD corrections 
to these processes are significant. We use the parton- level Monte Carlo program MCFM[p7|- 
ro| to estimate these backgrounds at NLO level. At the Tevatron, the probability P^jj that 
a jet fakes a photon is almost vanishing if the transverse momentum of the photon pJp is 
larger then 100 GeV because in this situation the hits in the central preradiator chambers 
are counted and the prompt photon is distinguished from meson decays|41]. However, at 



the LHC, to remain on the safe side, we set P^/j 



10 , as suggested in Ref.[^ 



Fig.|TT] shows the differential cross sections of both the signal and backgrounds as a func- 
tion of pTp and p^^^^. It can be seen that the Z'-f production is the main background. The 
distribution of the backgrounds decrease faster than that of the signal as the transverse 
momentum of the photon increases. Thus, the ratio of signal and background will increase 
if we set a larger px cut. 



Fig.|l2| shows the differential cross section of both the signal and backgrounds as a function 
of rj^. The main background is the Z7 production, and it is almost fiat in the full range of 
rj"'. In contrast, the signal concentrates in the central region of rj'^. This is a result of the 
scalar nature of the effective operator. These characteristics may help to select the events 
in experiments. 
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FIG. 11: Dependence of the differential cross section on p^ (left) and p™** (right). 
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FIG. 12: Dependence of the differential cross section on rj"' . 



Fig.O shows the integrated luminosity needed at the LHC with y/S = 14 TeV for 



5a discovery {S/VS + B = 5) of the signal. For A = 500 GeV and A = 1000 GeV, the 
integrated luminosities needed are around 0.03 fb~^ and 2.0 fb~^, respectively. The situation 
at the LHC with y/S = 7 TeV is also shown in Fig. 13. We find that the LHC may detect 
this signal once it collects an integrated luminosity of 1 fb~^, which means that we may 
expect the observation of this signal at the early stage of the LHC. From an experimental 
point of view, if we discover this signal then we can set an upper limit for the new physics 



scale, which has been illustrated in Fig. 14. 



16 




130 140 150 160 170 180 190 200 130 140 150 160 170 180 1«J0 200 



FIG. 13: The integrated luminosity needed at the LHC with \/5 = 14 TeV (left) and \/5 = 7 TeV 
(right) for a 5a discovery. The six curves from the bottom up correspond to the new physics scales 
of 500, 600, 700, 800, 900, 1000 GeV. 
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FIG. 14: The limits of the new physics scale for a 3cr exclusion and 5(T discovery at the LHC with 
\fS = 14 (7) TeV, assuming m = 150 GeV. 
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On the other hand, the LHC may not detect this signal at all. Thus we also present the 
exclusion limits of the new physics scale at the 3a [S /yB = 3) level in Fig.|T3. We can 



see that the new physics scale is constrained to be larger than 1450 (840) GeV if the LHC 
with vS* = 14 (7) TeV does not detect this signal after collecting an integrated luminosity 
of 10 ih-\ 

VI. CONCLUSION 

We have investigated DM annihilation and signal of DM and a photon associated pro- 
duction at the LHC induced by a dimension six effective operator at the NLO QCD level. 
We also study the main backgrounds from SM to this signal, i.e. Z boson and a photon 
associated production with invisible decay of Z boson, and Z boson and a jet production 
with the jet misidentified as a photon. We find that the p^ distributions of the backgrounds 
decrease faster than that of the signal with increasing of the transverse momentum of the 
photon. The rj'^ distributions of the backgrounds are almost flat in the full range of t]'^. 
In contrast, the signal lies mainly in the central region of t]'^. These characteristics may 
help to select the events in experiments. We show that in the parameter space allowed by 
the relic abundance constraint, which we have calculated at the NLO QCD level, the LHC 
with \/S = 7 TeV may discovery this signal at the 5a level after collecting an integrated 
luminosity of 1 fb~ . On the other hand, if this signal is not observed at the LHC, we set a 
lower limit on the new physics scale at the 3a level. 
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